arXiv:1504.06252vl [astro-ph.GA] 23 Apr 2015 


Mon. Not. R. Astron. Soc. 000 , 000-000 (0000) Printed 24 April 2015 (MN IATJ^X style file v2.2) 


Effect of turbulent velocity on the H I intensity 
fluctuation power spectrum from spiral galaxies 


Prasun Dutta 1 *, 

1 USER Bhopal, ITI Gas Rahat Building, Bhopal-62023, India. 


ABSTRACT 

We use numerical simulations to investigate effect of turbulent velocity on 
the power spectrum of H I intensity from external galaxies when (a) all emis¬ 
sion is considered, (b) emission with velocity range smaller than the turbulent 
velocity dispersion is considered. We found that for case (a) the intensity fluc¬ 
tuation depends directly only on the power spectrum of the column density, 
whereas for case (b) it depends only on the turbulent velocity fluctuation. We 
discuss the implications of this result in real observations of H I fluctuations. 
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1 INTRODUCTION 

Observationally, power spectrum of the H I intensity fluctuation in our Galaxy suggest ex¬ 


istence of scale invariant structures in 

the H I density over length sea 

es ranging as wide 

as sub parsec to a few 

rundred parsec ( 

Crovisier & Dickev 

1983; 

Green 

1993 

). These struc- 

tures are understood ( 

Elmegreen & Scalo 

2004 

) in terms of compressible fluid turbulence 


in the interstellar medium (ISM). In present theoretical understanding of ISM dynamics, 
compressible fluid turbulence plays an important role in the ISM evolution, energy trans¬ 
fer, star formation etc. Possible source of energy into the turbulence cascade is however 
debated, though is mostly ascribed to the supernova shocks as a large scale energy in¬ 
put. Different techniques have been developed to measure the velocity spectrum of the 
turbulence and hence infer the energy involved in the process. Techniques, originally devel¬ 
oped to get the velocity structures in the Galaxy, include statistics of centroid of velocities 
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reader here to VGA, which has been applied (IChepurnov et al. 


2009 


Chepurnov et al. 


2008 


(Es quivel fe Lazarianl 120091) . velocity coo rdinate spectrum or VCS (IPogosvan fc Lazarian 
20091 ). velocity channel anal ysis or VC A (ILazarian & Pogosvanl (120001 ) . henceforth LP00), 
spectral correlation function ( Padoan et all 2001) etc. We particular l y bring attention of our 


Spicker & Feitzinger 


2010 


Che p urno v & L azar ian 


1988) for observations in our Galaxy as 


well as nearby dwarf galaxies like Large and Small Magellanic Clouds. It was found that the 
velocity fluctuations also follow a power law. Interested reader may have a look at the article 
LP00 for a complete description of VCA, here we outline the basic principle behind this anal¬ 
ysis. VGA aims to extract the velocity power spectrum by comparing the power spectrum of 
intensity averaged over the entire velocity range of observation with the same averaged over 
a relatively small velocity range, smaller than the expected turbulence velocity dispersion. 
Differential rotation of our Galaxy allows us to have a direct mapping of the velocity values 
in the position-position-velocity data cube with the line of sight distance to the observing 
cloud. This in turn let us estimate the three dimensional power spectrum of the H I density 
fluctuations from the Galaxy. However, turbulence velocity fluctuations change the velocity 
to distance mapping and hence also modifies the intensity power spectrum. This is precisely 
what_VCA explores. 


Begum et al.1 (20Q6) have used a visibility based power spectrum estimator to measure 


the intensity fluctuation power spectrum of the nearby dwarf galaxy DDO 210. They infer 
that the density power spectrum has a slope of —2.75 over a length scale range of 80 to 
500 pc. They used the VCA with their position-position-velocity data cube and inferred 


an upper limit to the slope of the velocity power spectrum. 


Dutta et al. 


( 2008 . 2009a,b) 


;ernal dwar: 


Dutta et al 


and spiral galaxies and have estimated 
( 2013 ); Dutta fc Bharadwaj ( 2013 ) has 


has extended these study to several ex 
the density power spectrum. Recently,_ x _ 

estimated the power spectrum of the 18 spiral galaxies from THINGS 0 sample and found 
that the power spectrum of column density follow a power law over the length scales ranging 
400 pc to 16 kpc considering the entire sample. Slope of the power spectra for most of these 
galaxies was found to be within —1.5 to —1.8. Generating mechanism of these large scale 
structures are yet to be understood. Measuring statistics of the velocity fluctuations would 
help us understand the dynamical phenomena responsible for these structures. 

We note here two main difference between the position-position-velocity data cubes of 


1 THINGS: The H I Nearby Galaxy Survey ( Walte r et al 


200811 . 
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the H I emission observation in the Galaxy and external spiral galaxies. Line of sight to 
the observations for the H I emission from our Galaxy is mostly along the plane of the 
disk, while, the external galaxies are mostly for the face on galaxies and the line of sight is 
perpendicular to the disk. For observations in our Galaxy, different velocity slices of the data 
cube can be considered to be at different distances but at the same angular direction in the 
sky, whereas, for the external galaxies, different velocity slices of the data cube originates 
from different parts of the galaxy. This suggests that it would not be wise to directly use 
the results obtained in LPOO while inferring observations from external galaxies. 

As there exist no direct position to velocity mapping for the external spiral galaxies, an 
analytical investigation on the effect of the turbulence velocity is not straight forward, we 
refer to the numerical methods here. In this letter we perform numerical simulation to access 
how the intensity power spectrum is modified with the velocity fluctuations for the spiral 
galaxies. Section (2) gives a brief outline to our approach and the section (3) describes the 
numerical investigation we have performed. Results and discussions are discussed in section 
(4). We conclude in section (5). 


2 MODELLING H I EMISSION FROM SPIRAL GALAXY 

We adopt a coordinate system centred at the H I cloud in concern (or the external galaxy) 
with the line of sight direction aligned to the z axis, such that 

r = (x,y,z) = ( R,z ), (1) 


where R = (x,y) is a two dimensional vector in the sky plane. At small optical depth limit, 


the specific intensity of radiation (iDraind 20111 with rest frequency u 0 originated from a gas 
at r having temperature T is given as 


I{R,v) = J 0 J dzn HI (r)<j>(v), (2) 

where v = c(z/ 0 — u), v is the frequency of observation, J 0 = 3hl fl^ 21 and <f>(v) is the line shape 
function: 

(v - v z {r)) 2 
2cr 2 (r) 


<f>(v) = f o exp 


( 3 ) 
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k h T 


m HI 


is the 


Here v z {r) is the line of sight component of the velocity of the gas and cjr) = 
thermal velocity dispersion. In practice, the observed specihc intensity is always averaged 
over a velocity width Sv around v , hence 


-* 1 rv+Sv /2 

I obs (R,v,Sv) = — / dv'I(R, v'). 

OV Jv—5v/2 

Clearly, 


(4) 


_lim I obs (R, v, Sv) = I 0 N m (R ), (5) 

OV^OO 

where, Nhi(R) = f dznniif) is the column density. In practice, as the emission from the 
galaxy falls of to zero beyond a certain velocity, say ±Au, it is sufficient to carry the inte¬ 
gration in the above equation over the range —An to Av. Here we consider that the galaxy 
has no overall motion. 

Compressible fluid turbulence in the ISM of the galaxies induce scale invariant fluctua¬ 
tions in the density as well as velocity. Power spectrum of the column density fluctuation is 
given as 

Pn hi (K) = j dXe~ iS jl (N HI (R + X)N HI (R)), (6) 

where the averaging is performed over all possible values of R and all directions assuming 
homogeneity and isotropy in the random fluctuations. We define the power spectrum of the 
observed intensity fluctuation as 

P(K, Sv) = j dX e- iR *(I obs {R + X, v, 8v)I obs (R, v, Sv)) (7) 

Here we have assumed the homogeneity and isotropy of the intensity held and that the 
intensity power spectrum is independent of the centroid of the velocity v. Clearly, 


Jim P(K, Sv) oc P Nhi (K). 

5v^-Av 


( 8 ) 


This has been used ex tensively in literature to estimate t 


spectrum of our Galaxy 


2006 


Dutta et al. 


Crovisier & Dickey 


1983 


2009a) and spiral galaxies flDutta et al.l 2013 1). The power spectra is found 


re H I column c 


ensity power 


Green 1993), external dwarf 


Begum et al 


to follow power laws indicating turbulence to be operational, hence 


Pn hi (K) = A n „,K° (9) 

The line of sight component of velocity v z {r) has component from the systematic rotation 
of the galaxy v Q (r), as well as form the random motion of the cloud because of turbulence 


2 


fob : Boltzmann constant, m^ji '■ mass of hydrogen atom 
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v T {f), i.e, v z (r) = v n (r) + v T {f). Power spectrum of the turbulent velocity component is 
also expected to follow power law 

P v t(K) = A Vt K?. (10) 

LP00 has investigated the nature of P(K, Sv ) in detail in order to estimate the modification 
of the H I power spectrum by turbulence. They show that for observed H I gas in our galaxy, 
for a > —3, 

lim P(K, Sv) oc K a+ h/ 2 , (11) 

5v<ctt 

while as we expect, 

lim P(K, Sv) oc K a oc P Nhi {K). ( 12 ) 

Sv^-ctt 

Hence, by estimating the power spectra in two different limits above one can infer the velocity 
power spectrum slope. This method, usually known as velocity channel analysis, has been 
used to estimate the velocity fluctuation power spectrum of H I in our Galaxy and nearby 
dwarf galaxies. 

ft is important to realise that the particular direct linear mapping between v z and z was 
exploited in VCS, is rather different when we consider H I emission from an external galaxy. 
Considering tilted ring model, in later case, v ci (f) depends on the galacto-centric radius, 
position and inclination angles. Moreover, at a given v with Sv < Av, only a part of the 
galaxy’s disk is visible. In this letter we attempt to see how the turbulent velocity modifies 
the H I power spectrum for the external spiral galaxies and investigate if a similar procedure 
as VCS can be adopted to estimate the velocity fluctuation spectrum. 


2.1 Simplifications 


In order to simulate the H I emission from the external galaxies we adopt the following 
simplifications. Note that in this work we are not interested to simulate all aspect of the 
H I emission from the external galaxies, rather we are interested in investigating the modifica¬ 
tion in the power spectrum due to the turbulent velocity, which justifies these simplifications. 


• In case of a spiral galaxy the average H I profile Wfr) varies with the galacto-centric 
radius as well as in vertical direction. This leads to modification in the H I power spectrum 


as discussed in 


Duttaetah 


( 2009ap . Here, we consider W{r) to be independent of r. It is to 


be noted that this simplification also means that we are assuming the galaxy’s disk to be 
thick. We shall discuss the effect of this in the conclusions section. 
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Figure 1. Plot showing map of (a) Column density, (b) line of sight velocity, (c) observed specific intensity for Sv = ctt /2. 
Part (d) is the integrated line profile of the simulated galaxy with a = —2.5, /3 = —2.5, i = 10°. 


• Systematic rotation of the galaxy u Q (r), depends on the inclination and position angle 
as well as galacto-centric radius. To simplify matter we assume here that the position angle 
and inclination angle ( i ) of the galaxy do not change with galacto-centric radius and adopt 
a flat rotation curve with tangential velocity vq. In such a case, we can write 


v n (r) = 


v 0 sin(i)x 


lx 2 + y 2 cos 2 {i) 

ISM is known to be in pressure equilibrium (see 


(13) 


Wolfire et ah 


(2003j) and references 


therein) and more than one temperature gas coexists in it. This means, in principle, we 
need to consider different temperatures at different part of the galaxy and hence a varying 
cr(r'). This would give rise to an additional fluctuation in the observed specific intensity. Here 
we assume that the gas across the galaxy is at a constant temperature and we adopt the 
temperature that of the cold gas. 


3 SIMULATION 


We divide the simulation volume into NGrid 3 individual grids (cubes) which represents 
individual H I cloud with associated uhi and v z . As discussed in the previous section, we 
have kept a constant everywhere across the simulation volume. It can be shown that for a 
thick cube the power spectrum of the three dimensional density distribution has the same 
slope of that of the projected density, like the column density here (Dutta et al. (2009aj), 
LP 00). Hence, we generate Snni and v T such that they follow Gaussian distribution with 

















Effect of turbulent velocity on the H I intensity fluctuation power spectrum from spiral galaxies 7 


power spectrum of slope a and (3 respectively. iDutta fe Bharadwail (120131 1 h ave estimated 


the amplitude of the ff I fluctuations for six galaxies of the ThIINGS sample ( Walter et ah 


2008). They found that the amplitude of the column density fluctuations are approximately 
l/10 t/l that of the mean column density for the galaxies, ffence, here we consider 


n H i = n 0 [1 + fn HI 5n H i{rj ], 

with f nHI = 0.1 and Srinrij), the random component because of turbulence. 


(14) 


Tamburro et al 


(120091) has estimated the H I velocity dispersion for the galaxies in the THINGS from 
Moment-II maps and found that the turbulence velocity dispersion cr? vary in the range 
~ 5 to 20 km sec -1 . On the other hand, flattening velocity_of the rotation curve for the 


same galaxies has the range ~ 100 to 200 km sec 1 (jde Blok et al 


2008 b Here we adopt 


< j t = f v v o, with f v = 0.1. Considering the gas at a temperature of ~ 500 K, we adopt the 
thermal velocity dispersion to be a = frv o, where fr = 0.01. Note that the actual value of 
Vq is unimportant here. 

In order to see the effect of the density and velocity fluctuations in P(K, Sv), we consider 
a = (—1.5, —2.5), (3 = (—1.5, —2.5) and all combinations of these. In literature the velocity 
and density fluctuations because of turbulence are assumed to be uncorrelated. Here we 
consider two cases, either n^i and v T are completely uncorrelated or completely correlated. 
Galaxies with higher inclination angles can have scale mixing in the projected direction. 
Rotational velocity effects also manifest more at higher inclination angle because of the 
sin(i ) factor in eqn. (13). We choose the inclination angle to be i = 10°. 

Figure (1) shows different aspects of the simulated galaxy for a = —2.5, f3 = —2.5, 
i — 10° and for the case when rim and v T are uncorrelated. Column density map is shown 
in Figure (la), while the line of sight velocity v z is shown in Figure (lb). Figure (lc) shows 
I°hs f or a cer t a j n value of v with Sv = o' r /2. In this case only a part of the galaxy is visible 
and the area over which P(K, Sv) can be estimated is restricted. The integrated line profile 
of the galaxy is shown in Figure (Id). 

Assuming the centre of the galaxy to be at the centre of the simulation volume and 
the inclination angle to be i, we generated the specific intensity given in eqn. (1) with 
velocity resolution same as thermal velocity dispersion. Range of v is chosen such that all 
the emission from the model galaxy is included. We estimate P(K,8v ) defined in eqn. (7) 
for (a) Sv covering the entire H I emission, i.e, Sv = Av and (b) Sv = cqr/2. Results are 
discussed in the next section. 
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Figure 2. 2D power spectra with Sv = Av plotted for different combinations of a = (—1.5,—2.5) and = (—1.5, 2.5) for 
inclination angle i = 10°. Blue dash and dot-dash lines show the power spectra for /3 = —1.5 and —2.5 respectively with 
a = —2.5. The solid blue line is a plot of P(K) oc k a . The black curves are for corresponding power spectra with a = —1.5. 
All curves are shifted arbitrarily in the vertical direction for clarity. 


4 RESULTS AND DISCUSSIONS 

We first discuss the results for the case when we assume that the rim and v T are uncorre¬ 
lated. Figure (2) shows the power spectra P(K,Sv ) with Sv = Av. Here the “dot-dash” lines 
corresponds to the power spectra for j3 = —2.5 while the “dash-dash” lines are for f3 = —1.5. 
Power spectra corresponding to a = —2.5 are shown in blue and a = —1.5 are shown in 
black. The solid lines corresponds to power law with slopes —2.5 and —1.5 respectively. All 
curves are shifted in vertical direction arbitrarily for clarity. It is clear that the power spectra 
of the intensity with Sv = Av has the same slope that of the rim- irrespective of the slope 
of the velocity power spectra. As in our simulation we have considered a thick disk for the 
galaxy, we expect the slope of the power spectrum of column density to be same with that 
of uhi • Hence, nature of the power spectra in Figure (1) is quite expected and is just a 
verification of the eqn (8). 

We estimate the power spectra P(K, Sv) for all four combinations of a and /3 and i = 10° 
with Sv = cr-r/2, where we expect to see the effect of turbulence velocity v T in the intensity 
power spectra. Figure (3) shows the corresponding power spectra as in Figure (1) with 
Sv = a?/ 2. Before we interpret these curves, we need to realise that here emission is coming 
from only a part of the galaxy’s disk, as shown in Figure (lc). In such a case, the observed 
intensity power spectra would have effect of the shape of the window where the emission is 
coming from []. This is precisely why all four curves in Figure (3) have similar nature for 


3 Effect of the window is discussed in detail in 


Dutta et al. 


( 2009ah 
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Figure 3. 2D power spectra with Sv = crj-/2 plotted for different combinations as in Figure (1). Note that slope of the spectra 
changes for K > 0.8. 


K <0.8 and we can only expect to see the effect of a or ft beyond that. Interestingly, for 
K > 0.8 the power spectra is independent of the values of a and is different for different 
ft. This can be the effect of velocity modification, i.e, effect of the line of sight component 
of the turbulent velocity v T on the intensity power spectrum. As this is independent of a, 
the nature of the velocity modification is different than what is expected from the result of 
LP00 (see eqn. (10, 11)). 

To investigate how P(K, Sv) for K > 0.8 changes with different values of ft, we performed 
the same simulation with a = —2.5 and i = 10°, for values of ft ranging —3.0 to —1.0. For 
each case, we fit the power spectra at K > 0.8 with a power law of the form P(K) oc K 1 and 
note the best fit values. Since in simulation we have not added any contribution from the 
observational uncertainties, we only use the sample variance generated noise to do this fit. 
Results are shown in Figure (4), where we have plotted ft in x axis with 7 with errors form 
the fit in y axis. We use a second order polynomial to empirically fit the values of 7 against 
(3, i.e, f(x) = a 0 + aix + a 2 x 2 with a 0 ,a\ and a 2 values —0.27, —0.17 and —0.15 respectively. 

Next we consider the case when the fluctuations in rim and that in v T are correlated, 
in fact we use the same set of gaussian random variables to represent them. Hence, in this 
case, we only consider variation of a since nni and v T fluctuations are scaled up version of 
the same original gaussian variables. As expected the power spectrum with Sv = Av has the 
slope of the column density power spectrum and the corresponding plot is exactly similar 
to Figure (2) except from minute differences arising because of statistical fluctuations. We 
estimated the power spectrum of H I intensity with Sv = a T / 2. These power spectra also 
show trends as in Figure (3), for K < 0.8 it is dominated by the windowing effect and for 
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Figure 4. Points with error bars represents best fit values of 7 as a function of ft. We use a polynomial fit to the points and 
the best fit line is shown in blue solid line. The values of the polynomial coefficients are also given in the plot window. 


larger K, it is a power law with exactly similar variation of the slope of the spectra as in 
Figure(4). We do not show these plots here to avoid repetition. 

To summarise, for both the cases uhi and v T uncorrelated and perfectly correlated the 
power spectrum P(K,8v ) with Sv = Av, we always reproduce power law with the same 
slope as n H i power spectrum and with Sv = <tt/2, at larger K the power spectrum has a 
certain slope that only correlates with the slope of the velocity spectrum. 


5 CONCLUSIONS 


In this letter we investigate how the H I intensity fluctuation power spectrum is related to 
the number density and the line of sight component of the turbulent velocity for external 
galaxies. We found that for scale invariant fluctuations in both density and velocity, when 
the emission is integrated over all velocity range, the intensity fluctuation power spectra 
follow a power law that has the same slope of the power spectra of the H I number density 
fluctuations. We consider a thick disk for the galaxy in this case, in case of thin disk, the 
intensity spectra would have slope shallower by order unity f lDutta et al.l (2 009a ). LP00). 

When the emission is integrated over a velocity range smaller than the turbulent velocity 
dispersion, due to the galactic rotation, only a part of the galaxy is visible. Effect of the 
density or the velocity fluctuation in the intensity power spectra can be inferred only for 
higher values of K and for a relatively narrow range of K values. We found that the slope of 
the spectra at these range follow approximately a power law with slope 7 related nonlinearly 
only to the slope of the velocity power spectrum /3 and independent of the power spectrum of 
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the density. This is a different result compared to LPOO, where it is expected that 7 = a+f3/ 2. 
Note that, this result is based on a power law fit to the power spectrum for a shallower 
range of K. Nevertheless, it clearly demonstrates that velocity modification of the H I power 
spectrum for external galaxies is quite different from that in our Galaxy. 


In our simulations, we did several simplifications. 


of the galaxy. However, as it is shown in 


Dutta et al. 


irst we ignored the overall H I profile 


( 2009al ). effect of this profile is to mod¬ 


ify the power spectra at lower K , which considering the galaxy spanning over the simulation 
volume would happen at K <0.1 and is not the regime of interest for the velocity modi¬ 
fication. Similarly, including a realistic rotation curve would have only change the window 
over which the power spectrum can be investigated for velocity modification. Given these, 
our results also stand out for real galaxies. Effect of varying thermal velocity dispersion over 
the galaxy, on the other hand, is more complicated and needs to be investigated in detail 
separately. 


Finally, we discuss the feasibility to use the relation we obtain between 7 and f3 in 
Figure (4) for a real observation. Considering the galaxy spread over our entire simulation 
volume, the dynamic range in K from simulation is approximately same as that in the 


THINGS observation. 


Dutta et al. 


(120131) has estimated the power spectra of 18 nearby 


spiral galaxies from THIGNS sample. Given the baseline coverage of observation, they could 
estimate the power spectra till l/4 th the largest baseline. This is because at higher baselines 
the baseline coverage is restricted and signal to noise is insufficient to estimate the power 
spectrum with statistical significance. We use NGrid = 512 for performing our simulation. 
As the largest K in simulation is 7 r, l/4 th the available baseline rang of THINGS corresponds 
to K = 7 t/4 ~ 0.8 in our simulation, making these observation inefficient to probe the 
velocity modification this way. We choose an inclination angle of 10° for our simulation, any 
higher inclination angle would result rather even smaller range in K over which the values 
of 7 can be estimated. 


We conclude that with present telescopes to use VGA techniques for external galaxies, we 
need high integration time. An alternate method to estimate the turbulent velocity spectra 
of the external galaxies would be more useful. We aim to investigate along this direction in 
future. 
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